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基于 自 适应 吸引 半径 的 萤火虫 算法 的 粒子 滤波 


EPE la, lb. 潘 Ed la, 1b 
d. 首都 师范 大 学 a 信息 工程 学 院 ; b. 北京 成 像 技术 高 精 炎 中 心 ， 北京 100048) 


摘 ， 要 : 针对 粒子 滤波 算法 对 粒子 数目 的 大 量 需 求 等 嗓 端 ， 提 出 一 种 基于 改进 的 莹 火 虫 算法 的 粒子 滤波 。 首 先 ， 在 莹 
火 吕 的 亮度 公式 中 引入 观测 值 信 息 ， 以 提高 算法 跟踪 的 准确 性 ; 其 次 ， 提 出 自 适 应 吸引 半径 参数 来 控制 营 火 忠 群 寻 优 
时 的 吸引 范围 ， 使 算法 的 实时 性 更 好 ; 最 终 利用 萤火虫 算法 的 选 代 寻 优 来 进行 粒子 更 新 。 对 比 实验 表明 ， 该 算法 在 跟 
踪 精 度 和 运行 时 间 上 都 有 所 优化 ， 说 明 该 算法 即使 在 粒子 数 目 较 少 的 条 件 下 ， 也 能 保证 目标 跟踪 的 准确 性 和 实时 性 。 
关键 词 ， 自 适应 吸引 半径 ; ARE: 目标 跟踪 ; 粒子 多 样 性 :相对 亮度 
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Particle filter based on firefly algorithm with adaptive attraction radius 
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University, Beijing 100048, China) 


Abstract: For the large demand of particles when particle filter tracking, this paper proposed a particle filter algorithm based 
on the improved firefly algorithm. Firstly, in order to improve the accuracy of tracking, the improved algorithm introduced the 
observation value into the brightness formula of the firefly; Secondly, when the fireflies swarm optimized, an adaptive attraction 
radius parameter was proposed to control the attraction range, so that the real-time performance would be better; Finally, the 
iterative optimization ofthe firefly algorithm was used to update the particles. Comparison experiments show that, the algorithm 


has been optimized for tracking accuracy and running time. It shows that the algorithm can guarantee the accuracy and real-time 


performance of the target tracking even if the number of particles is small. 
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致 算法 实时 性 差 ， 重 采样 后 造成 的 粒子 贫 化 问题 。 因 此 ， 基 于 
群 智能 优化 算法 的 粒子 滤波 是 一 个 轿 新 的 发 展 方向 。 这 类 算法 
人 们 在 生活 中 所 获取 到 的 信息 超过 80% 来 源 于 视觉 ， 随 着 。 通过 对 粒子 分 布 进行 迭代 寻 优 减少 对 粒子 数 的 需求 ， 同 时 算法 


0 引言 


计算 机 视觉 领域 的 飞速 发 展 ， 基 于 视频 序列 的 目标 跟踪 作为 计 不 涉及 对 低 权 值 粒子 的 舍弃 ， 从 根本 上 解决 了 粒子 贫 化 问题 。 
算 机 视觉 应 用 中 的 重要 任务 ， 被 广泛 应 用 在 视频 监控 、 车 辆 跟 张 才 千 等 人 由利 用 粒子 群 算法 (PSO) 在 重 采样 前 对 粒子 集 
踪 、 手 势 识 别 等 众多 领域 0”1。 进行 优化 ， 使 粒子 尽 可 能 多 地 分 布 在 高 似 然 区 域 ， 有 效 地 防止 


传统 的 目标 跟踪 主要 是 采用 基于 模型 线性 化 和 高 斯 假设 的 了 重 采 样 后 粒子 多 样 性 的 降低 。Wang 和 Heris 等 人 M9 分 别 在 
卡尔 曼 滤 波 来 完成 ， 但 在 实际 应 用 中 ， 会 存在 大 量 的 非 线性 系 ” 重 采 样 阶段 和 初始 化 采样 阶段 使 用 蚁 群 算法 来 优化 粒子 滤波 进 
统 或 非 高 斯 噪声 系统 ， 卡 尔 曼 滤 波 显 然 已 经 满足 不 了 系统 要 求 行 目标 跟踪 ， 虽然 很 好 地 解决 了 粒子 匮乏 及 粒子 退化 的 问题 ， 
该 方法 并 没 


B]。 基 于 蒙特 卡 罗 方 法 的 粒子 滤波 外 由 于 在 跟踪 过 程 中 不 需要 但 由 于 搜索 过 程 和 结束 条 件 带 有 很 强 的 经 验 性 ， 有 

做 线性 高 斯 的 假设 ， 可 以 很 好 地 求解 非 线性 非 高 斯 系统 下 的 跟 ” 得 到 很 好 的 使 用 。Wang 等 人 0 使 用 改进 的 遗传 算法 作为 重 采 

踩 问题 ， 因 此 得 到 了 广泛 的 应 用 5 T. 样 方法 ， 设 计 交 叉 变 异 概率 随 粒 子 的 退化 程度 成 正比 变化 来 控 
然而 ， 尽 管 粒子 滤波 在 目标 跟踪 方面 有 其 一 定 的 优越 性 ， 制 交 叉 变异 操作 采样 新 的 粒子 ， 虽 然 提 高 了 粒子 的 有 效 性 和 多 

例如 可 以 利用 蒙特 卡 罗 方 法 将 积分 运算 转换 为 求 和 运算 等 , 但 ， 样 性 ， 但 其 迭代 次 数 和 浆 值 不 易 确 定 。 

它 存在 着 一 些 缺 点 ， 如 需要 大 量 的 粒子 来 保证 跟踪 的 准确 性 导 相 比 于 以 上 几 种 算法 ， 莹 火 虫 算法 作为 最 新 的 群 短 能 优化 
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算法 之 一 ， 具 有 较 高 的 收敛 精度 、 


更 快 的 收敛 速度 且 容 易 实现 


的 优点 ， 


但 目前 国内 外 关于 萤火虫 算法 优化 粒子 滤波 的 工作 还 


民 少 。 朱 文 超 等 人 0 根据 权 值 的 退化 程度 将 粒子 集 分 层 并 将 优 
化 层 粒子 映射 到 高 似 然 区 域 ， 但 不 能 真正 解决 粒子 贫 化 问题 。 


Gao 等 人 [3 实现 了 改进 的 吸引 度 公 


S 式 的 萤火虫 算法 ， 使 粒子 集 


分 布 到 高 似 然 区 域 来 优化 粒子 滤波 ， 但 亮度 计算 公式 的 特性 会 


导致 很 高 的 运算 复杂 度 
子 滤波 (FA-PF)， 使 ) 
寻 优 过 程 中 出 现 局 部 极 值 的 现象 ， 同 时 也 有 效 地 控 
子 多 样 性 的 衰减 ， 但 


| 梦 楚 等 人 0 提出 莹 火 虫 算法 优化 粒 
] 改 进 的 位 置 更 新 公式 进行 迁 代 寻 优 ， 避 旬 
制 了 粒 
于 粒子 相互 吸引 时 需要 全 部 粒子 之 间 两 


两 进行 交互 运算 ， 所 以 会 造成 计算 复杂 度 的 增加 和 运行 时 间 的 


针对 上 述 问题 
火 虫 算法 的 粒子 滤波 ， 通 
子 相 互 吸引 时 的 吸引 范围 ， 


， 本 文 提出 了 一 种 基于 自 适 应 吸引 半径 的 萤 
过 加 入 自 适 应 吸引 半径 参数 来 控制 粒 
解决 了 粒子 交互 时 参与 运算 的 粒子 


N 


数目 较 大 以 及 计算 复杂 度 增加 的 问题 。 
1 ”粒子 滤波 算法 


粒子 滤波 算法 
噪声 的 观测 值 来 估计 系统 状态 的 后 验 密度 
系统 状态 的 估计 09。 粒 子 的 运动 状态 是 一 阶 自 回 归 状 态 
状态 模型 和 观测 模型 数学 表示 如 下 : 


zi, 


其 中 : fA h 分别 为 状态 函数 和 观测 函数 ; 
统 噪声 和 观测 噪声 。 


所 示 。 


,sG hinaXiv 合 作 期 天 


算法 的 粒子 滤波 


新 计算 如 下 : 


Py, | X * Yua P(x, | Yua) 


i | »4)- P(y, | Yia) P (Yua) 


(4) 
P(y, | x, )P(x, | Jua) 
P(y, | Yaa) 


以 上 计算 可 知 ， 在 进行 粒子 预测 和 
到 所 有 的 观测 数据 ， 随 着 时 间 的 推移 ， 会 造成 计算 量 的 大 大 


更 新 的 过 程 中 ， 需 要 


增加 。 粒 子 滤波 采用 序 货 重 要 性 采样 (SIS) 的 方式 进行 递 轨 计算 ， 
Tidi k —1 时 刻 的 粒子 权 值 计算 得 到 k 时 刻 每 个 粒子 的 重要 性 
权 值 如 式 (5) 所 示 ， 并 进行 权 值 归 一 化 如 式 (6) 所 示 。 
kay pP(ydx) px) 
"i E q(x Xok-1> yx) 6) 
Or 
Yw (6) 


izl 


基于 上 述 内 容 , 粒子 滤波 在 k=0 时刻 对 所 有 粒子 进行 初始 
化 状态 信息 ， 随 着 时 间 的 推移 不 断 地 对 粒子 的 状态 进行 预测 和 
更 新 ， 计 算 所 有 粒子 的 权 值 并 归 一 化 ， 判 断 是 否 需 要 重 采样 。 


的 关键 是 使 用 一 组 带 权 值 的 粒子 集 以 及 带 有 
函数 53， 从 而 实现 对 


空间 模 


hf (md) 0) 
», =h(mv) Q) 
ux 和 vi 分 别 为 系 
状态 模型 和 观测 模型 的 动态 表示 0 如 图 1 


图 1 状态 模型 和 观测 模型 的 动态 系统 


在 已 知 k-1 时 刻 的 粒子 状态 条 件 下 ， 根 据 状态 模型 、 观 测 


模型 以 及 贝 叶 斯 估计 理论 对 大 时 刻 的 粒子 状态 进行 预测 和 更 新 
09。 粒 子 状态 预测 的 计算 如 下 : 


但 


积分 运算 转换 为 有 限 样 


刻 的 观测 


P(x | x Yu Pa wa) 
P(yua) 


dx 


k-l 


P(x | Yua) =Í 
= f PG, | xi Yua )P(x,, , Yia) dx, , 6) 


= [PG 1) P Qua 1) d 
于 积分 运算 不 易 实 现 ， 粒 子 滤 波 采用 蒙特 卡 罗 方 法 将 
本 点 的 求 和 和 运算。 预测 之 后 ， 根 据 大 时 
值 及 贝 叶 斯 估计 对 粒子 状态 进行 更 新 。 粒 子 的 状态 更 


若 需 要 ， 则 进行 重 采 样 操作 (复制 权重 大 的 粒子 ， 舍 弃权 重 小 
的 粒子 ) 得 到 一 组 新 的 粒子 的 状态 ， 若 不 需要 ， 则 进入 下 一 时 
刻 的 预测 更 新 ， 实 现 对 目标 的 循环 跟踪 。 
2 ”萤火虫 算法 

群 智能 优化 算法 是 通过 对 生物 种 群 的 群 行为 进行 数学 建 模 
来 达到 寻 优 目的 的 一 类 算法 。 董 火 虫 算法 作为 群 智 能 优化 算法 
的 一 种 ， 主 要 是 模仿 自然 界 中 萤火虫 的 发 光 行 为 ， 以 个 体 间 的 
亮度 差异 为 基础 进行 相互 吸引 完成 寻 优 过 程 ， 达 到 寻 优 目的 。 
pup 主要 依靠 亮度 和 吸引 度 


萤火虫 之 间 的 相互 吸引 
21 萤火虫 的 相对 亮度 

萤火虫 的 相对 亮度 是 萤火虫 之 间 相 互 吸引 的 基础 。 在 相互 
吸引 的 过 程 中 ， 由 于 萤火虫 两 两 个 体 的 亮度 不 同 ， 亮 度 较 高 莹 
火 虫 会 吸引 亮度 较 低 萤火虫 使 莹 火 虫 种 群 一 致 向 最 亮 莹 火 虫 移 
动 来 完成 寻 优 过 程 。 相 对 亮度 公式 如 下 : 

I =e" (7) 

其 中 : 鳌 为 草 火 虫 的 相对 亮度 大 小 ; Ho 为 最 亮 草 火 虫 的 亮度 ， 
最 优 董 火 虫 所 处 位 置 的 目标 函数 值得 到 ; 7 为 光 强 吸收 系数 ， 
主要 用 来 体现 亮度 随 距 离 的 增加 而 减弱 的 特性 ， 一 般 取 
y e[0.01,100] ， 方 为 萤火虫 ?到 萤火虫 7 的 篆 卡 儿 距 离 。 
2.2 萤火虫 的 吸引 度 

在 萤火虫 相互 吸引 的 过 程 中 ， 通 过 吸引 度 来 衡量 
体 吸 引力 的 大 小 ， 结 合 其 个 体 的 相对 亮度 ， 
引 度 公式 如 下 


萤火虫 个 
完成 位 置 更 新 。 吸 


B, 7 Ae (8) 
其 中 ， ,为 最 大 吸引 度 ， 即 在 光源 处 莹 火 虫 的 吸引 度 。 
2.3 ”萤火虫 的 位 置 更 新 
萤火虫 种 群 中 的 个 体会 在 吸引 与 被 吸引 的 过 程 更 新 自身 的 
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C 


LES E. Ma AS BHL TA re K R E EPRE fa Eo RES 


高 的 区 域 。 位 置 更 新 公式 如 下 : 


LET * f(x, 一 总 )+ ae 


(9) 
腿 设 萤火虫 了 被 萤火虫 ;吸引 ， 莹 火 虫 了 将 按 如 上 公式 进 
TME EH. Feb, xox; HEKE i j 的 位 置 ，B 为 萤火虫 i 对 


萤火虫 j 的 吸引 度 ; e 为 常数 ， 一 般 取 Qe[0,1; 6; 是 由 高 } 


分 布 、 均 匀 分 布 或 者 其 他 分 布 得 到 的 随机 数 。 
3 ”基于 改进 的 萤火虫 算法 的 粒子 滤波 


3.1 算法 描述 
本 文 主要 从 以 下 两 个 方面 对 算法 进行 了 改进 : 


及 粒子 间距 离 来 计算 相对 亮度 ， 若 将 该 相对 亮度 直接 应 用 于 与 
粒子 滤波 结合 完成 粒子 的 迭代 寻 优 过 程 ， 无 法 实时 地 利用 粒子 
最 新 的 观测 值 信息 ， 导 致 算法 的 准确 性 不 高 。 同 时 ， 其 计算 方 
式 的 特殊 性 会 在 算法 运行 时 带 来 巨大 的 计算 量 。 

针对 该 缺陷 ， 本 文采 用 融入 粒子 最 新 的 观测 值 信息 的 方法 
来 提高 算法 准确 性 ， 减 少 运行 时 的 计算 量 。 改 进 后 的 相对 亮度 
数学 表达 如 下 : 


1=(y.—») (10) 
其 中 : 了 为 改进 后 的 相对 亮度 Xi 为 粒子 的 观测 值 ，3 为 粒子 
的 预测 值 。 若 将 相对 亮度 公式 扩充 到 多 维 情况 中 ， 关 键 是 计算 
观测 值 和 预测 值 在 多 维 空间 中 空间 差 值 的 平方 。 对 于 二 维 坐标 


a) 由 于 原始 的 亮度 公式 计算 复杂 度 较 高 ， 并 且 没 有 结合 粒 


子 最 新 的 观测 值 信息 
本 文 将 观测 值 信息 引入 到 萤火虫 算法 
的 准确 性 。 


T 


\, 这 会 导致 跟踪 准确 性 下 降 的 问题 ,因此 ， 
的 亮度 公式 中 来 提高 跟踪 


bp) 考虑 到 粒子 数目 增加 时 ， 粒 子 集 在 进行 相互 吸引 寻 优 的 


过 程 中 ， 计 算 量 将 呈 指 数 增长 的 问题 。 本 文 在 充分 研究 董 火 虫 


信息 ， 相 对 亮度 公式 应 为 
I2(y,- y) +- (11) 
其 中 : Guy 是 目标 观测 坐标 值 ，(x,y) 是 目标 预测 坐标 值 。 
相对 于 原始 的 相对 亮度 公式 而 言 ， 改 进 后 的 公式 利用 每 个 
时 刻 每 个 粒子 自身 的 观测 值 与 预测 值 之 间 的 差 值 平方 来 表示 粒 
子 每 个 时 刻 的 相对 亮度 ， 不 再 需要 将 全 部 粒子 与 最 优 粒 子 进行 


算法 特点 的 基础 上 引入 一 个 新 的 参数 一 吸引 半径 。 它 主要 是 随 
粒子 的 亮度 变化 ,可 以 有 效 地 缩减 粒子 交互 时 的 粒子 吸引 范围 。 
与 文献 [13] 相 比 ， 本 文采 用 融入 观测 值 信息 的 相对 亮度 作 


为 萤火虫 算法 迁 代 寻 优 的 基础 ， 使 得 用 于 跟踪 目标 的 粒子 集 能 
够 更 加 集中 地 分 布 在 高 似 然 函数 区 域 ， 跟 踪 更 加 准确 。 同 时 当 
多 数 粒子 具有 高 权 值 时 ， 表 示 对 准确 跟踪 有 贡献 的 粒子 占据 了 


粒子 集 的 绝 大 部 分 ， 而 对 跟踪 无 明显 意义 的 小 权 值 粒子 数目 


大 减少 ， 造 成 参与 重 采样 步骤 的 粒子 数量 有 了 大 幅度 的 缩减 ， 


比较 计算 ， 减 少 了 算法 在 欠 代 寻 优 过 程 中 的 计算 量 。 同 时 ， 通 
过 加 入 粒子 观测 值 信息 的 策略 ， 将 算法 预测 值 与 观测 值 的 差 值 
实时 地 反映 在 相对 亮度 中 ， 而 相对 亮度 的 变化 会 影响 下 一 时 为 
粒子 的 预测 值 的 估计 ， 从 而 有 效 地 提高 了 算法 跟踪 的 准确 性 。 
3.2.2 自 适 应 吸引 半径 的 提出 

萤火虫 算法 的 寻 优 过 程 是 对 全 部 草 火 虫 个 体 的 亮度 进行 两 
两 相互 比较 ， 亮 度 较 低 的 萤火虫 会 被 亮度 较 高 的 莉 火 虫 吸引 ， 
更 新 自身 的 位 置信 息 和 亮度 信息 。 但 若 以 原始 吸引 方式 进行 


从 而 减少 了 算法 在 重 采 样 过程 中 的 时 间 损 耗 ， 有 效 地 缩减 了 算 


法 跟踪 时 的 运行 时 间 。 


本 文 算法 在 进行 目标 跟踪 过 程 中 的 粒子 描述 如 图 2 所 示 。 


o oo 初始 化 粒子 


预测 


萤火虫 算法 


重 采 样 


图 2 ”基于 萤火虫 算法 的 粒子 滤波 的 粒子 描述 区 
3.2 ”基于 改进 的 萤火虫 算法 的 粒子 滤波 


针对 粒子 滤波 跟踪 时 所 需 的 粒子 数目 庞大 的 问题 ， 基 于 群 


寻 优 ， 随 着 萤火虫 群 数目 的 增加 ， 会 造成 计算 量 的 巨 幅 增长 ， 
带 来 很 大 的 时 间 损 耗 。 因 此 ， 本 文 设 计 一 个 新 的 参数 一 一 吸引 
半径 。 吸 引 半 径 的 作用 主要 是 限制 萤火虫 个 体 间 相互 吸引 时 的 
吸引 范围 以 及 粒子 数目 。 吸 引 半 径 的 大 小 由 亮度 决定 ， 亮 度 越 
大 的 粒子 吸引 半径 越 大 。 

萤火虫 亮度 与 吸引 半径 的 关系 如 图 3 所 示 。 


D 
a 


cV 


43 ”亮度 与 吸引 半径 关系 
a,b,c 分 别 表示 三 个 亮度 不 同 的 萤火虫 ， 亮度 为 4a>b>c， 


智能 优化 算法 的 粒子 滤波 是 目前 解决 该 问题 的 主要 方法 。 但 


于 粒子 群 算法 需要 存储 的 数据 量 大 ， 不 易 更 新 ， 而 蚁 群 算法 


的 


寻 优 过 程 完全 依赖 于 过 程 中 反馈 的 信息 ， 速 度 慢 ， 所 以 本 文 提 


出 基于 自 适 应 吸引 半径 的 萤火虫 算法 的 粒子 滤波 。 
3.2.1 萤火虫 相对 亮度 公式 的 改进 
原始 的 萤火虫 相对 亮度 公式 即 式 (7) 仅 依靠 最 亮 粒子 信 


显而易见 , 吸引 半径 为 六 > 态 > 天 。 对 于 萤火虫 ,由 于 六 > 各 
RIS A HR b des Jed a 的 吸引 范围 内 ， 所 以 萤火虫 会 被 萤 火 
虫 & 吸引 并 更 新 位 置 ， 但 对 于 萤火虫 < ， 尽 管 其 亮度 最 小 ， 但 
于 它 不 在 萤火虫 % 的 吸引 范围 内 , 便 不 会 被 萤火虫 a,b 吸引 
更 新 位 置 。 因 此 ， 对 于 两 个 亮度 不 同 的 萤火虫 ， 只 有 当 莉 火 
虫 间距 离 小 于 较 亮 草 火 虫 的 吸引 半径 时 才 会 发 生 位 置 更 新 ， 否 
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录用 定稿 
则 不 更 新 位 置 。 

在 本 文中 采用 粒子 观测 值 与 预测 值 差 值 的 平方 来 表示 粒子 
的 相对 亮度 ， 则 相对 亮度 值 越 小 ， 表 示 粒 子 的 预测 值 和 观测 值 
越 接 近 ， 越 靠近 最 优 粒子 ， 粒 子 的 吸引 半径 则 应 该 越 大 。 设 计 
吸引 半径 计算 公式 为 

r=1% (12) 

其 中 ; 为 粒子 的 吸引 半径 ， 了 为 粒子 的 亮度 大 小 。 


基于 自 适应 吸引 半径 的 
时 ， 不 再 需要 将 每 一 个 粒子 与 其 他 全 部 粒 了 


莹 火 虫 算法 在 进行 粒子 间 相 互 吸引 
进行 亮度 比较 并 更 


新 位 置 ， 只 需 将 粒子 与 在 其 吸引 半径 内 的 粒子 进行 亮度 比较 完 
成 吸引 和 位 置 更 新 过 程 即 可 。 这 使 得 算法 在 迭代 寻 优 过 程 中 参 


与 计算 的 粒子 数目 大 大 降低 ， 从 而 减少 了 算法 准确 跟踪 所 需 的 
粒子 数目 ， 有 效 地 降低 了 算法 的 计算 复杂 度 ， 缩 减 了 算法 的 运 
行 时 间 。 

3.3 ”算法 步骤 


算法 1 基于 自 适 应 吸引 半径 的 萤火虫 算法 的 粒子 滤波 
k =0 时刻 初始 化 
if Kk «Kk. 
fo i<N 
fo j«N 


计算 天 = 大 二 1 时 刻 粒子 预测 值 
计算 粒子 i 的 相对 亮度 及 粒子 i 的 吸引 半 
计算 粒子 j 之 间 的 相互 距离 广 


长 
HN 
E 


if rur 
计算 粒子 J 的 位 置 更 新 
end 
j=j+1 
end 
i=i+1 
end 
end 
更 新 权 值 并 归 一 化 
if Ng SN 
end 
状态 估计 
结束 
传统 粒子 滤波 算法 使 用 贝 叶 斯 估计 理论 对 粒子 集 进行 更 新 
本 文 算法 引入 了 基于 自 适应 吸引 半径 的 萤火虫 算法 更 新 粒子 集 
信息 。 在 获得 粒子 的 预测 值 和 观测 值 信息 后 ， 利 用 萤火虫 算法 
对 粒子 集 进行 迭代 寻 优 并 完成 粒子 状态 的 更 新 ， 使 粒子 集 分 布 


在 高 似 然 区 域 ， 
4 ”实验 与 分 析 


实验 条 件 
实验 软 硬 件 环境 见 表 1。 


进行 更 准确 的 跟踪 。 


4.1 
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算法 的 粒子 滤波 


表 1 实验 软 硬 件 环 境 
硬件 环境 软件 环境 
Intel? is 处 理 器 MATLAB R2016a 
4.00 GB 内 存 Windows 10 操作 系统 


实验 采用 的 系统 状态 模型 和 观测 模型 如 下 : 


sosete ese) co) (13) 
»()-20 *v(r) (14) 


该 系统 为 单 变量 非 静态 增长 模型 ， 是 一 个 典型 的 非 线性 系 
统 ， 传 统 的 粒子 滤波 方法 很 难 进行 准确 的 跟踪 。 设 观测 噪声 方 
差 R=1， 滤 波 时 间 长 度 为 50。 在 萤火虫 算法 中 ， 一 般 最 大 吸引 
度 的 设置 区 间 为 [0.8,1]， 本 文 设置 最 大 吸引 度 为 0.8， 光 强 吸收 


系数 设置 为 1。 
42 定量 分 析 

在 软 硬 件 条 件 相 同 的 情况 下 ， 对 比 算法 使 用 标准 的 粒子 滤 
波 算法 (以 下 简称 PF) 和 文献 [14] 中 提出 的 萤火虫 算法 智能 优化 
粒子 滤波 算法 (以 下 简称 FA-PF)， 记 录 不 同 算法 的 运行 结果 并 


进行 比较 和 分 析 。 


本 文 所 有 结果 数据 都 是 通过 100 次 重复 实验 取 平 均 数 据 值 
得 到 的 。 
4.2.1 精度 测试 

本 文 使 用 均 方 根 误差 (又 称 标准 误差 , RMSE) 公 式 作 为 误差 


大 小 即 预测 精度 的 衡量 方式 。 均 方 根 误差 是 用 来 衡量 观测 值 与 
预测 值 之 间 的 偏差 ， 能 够 很 好 地 反映 出 预测 的 精密 度 。 均 方 根 
误差 的 计算 公式 为 


T 
RMSE = SIC) 


n 


(15) 


其 中 :了 为 跟踪 的 时 间 长 度 ， 名 为 状态 均值 。 

根据 均 方 根 误差 的 计算 公式 可 知 ， 状 态 估计 值 越 接近 于 观 
测 值 ，RMSE 的 值 越 小 ， 跟 踪 准 确 性 较 高 ， 算 法 性 能 越 好 ， 反 
之 则 越 差 。 

D 算法 性 能 测试 

在 不 同 粒子 数 


N=50、N=100、N=200 及 不 同 过 程 噪声 


Q-1. Q-0.1 条 件 下 运行 本 文 算法 与 PF、FA-PF 三 种 算法 ， 
记录 不 同 条 件 下 跟踪 的 均 方 根 误差 值 ， 结 果 数 据 见 表 2。 
K2 不 同 粒子 数目 下 几 种 算法 的 RMSE XJ EE 
参数 PF FA-PF ios uiid 
法 对 比 PF ”对 比 FA-PF 
N=50 4.9413 4.0585 3.8417 22394 5.396 
=]  N-100 4.3735 3.5640 3.327] 23.9% 6.696 
N-200 4.0221 3.4576 3.1543 21.694 8.896 
N-50 33104 32162 3.1439 5.1% 2.394 
Q-0.1 N-100 3.2053 3.1507 2.9825 6.9% 5.396 
N-200 3.0060 2.96903 2.8555 5.0% 3.896 


表 2 可 知 ， 本 文 算法 在 不 同 粒子 数 


、 过 程 噪声 条 件 下 
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RMSE 值 都 较 其 他 两 种 算法 小 ， 说 明 本 文 算法 的 准确 性 在 不 同 
情况 下 均 高 于 PF 和 FA-PF 算法 。 在 不 同 过 程 噪声 下 ， 随 着 粒 
子 数 目 N 的 增加 ， 三 种 算法 的 准确 性 都 有 所 提高 。 

除 此 之 外 ，Q=1 时 ， 本 文 算法 的 准确 性 在 N=100 时 依然 高 
于 FA-PF 算法 在 N=200 时 ; Q=0.1 时 , 本 文 算法 在 N=100 时 的 
准确 性 近似 于 FA-PF 算法 在 N=200 时 的 准确 性 , 说 明 本 文 算法 
在 粒子 数目 少 的 情况 下 依旧 能 保证 跟踪 的 准确 性 。 

究 其 原因 ， 由 于 本 文 算法 在 重 采 样 前 利用 改进 后 的 萤火虫 
算法 对 粒子 集 进行 迭代 寻 优 操作 ， 使 粒子 分 布 在 高 似 然 区 域 粒 
子 的 权 值 都 较 高 ， 进 行 准确 跟踪 所 需 的 粒子 数目 也 有 所 下 降 。 

通过 表 2 数据 计算 可 得 ， 在 不 同 的 过 程 噪声 下 ， 本 文 算法 
相对 于 PF 以 及 FA-PF 的 平均 提高 率 如 表 3 所 示 。 

表 3 不同 过 程 噪声 下 本 文 算法 准确 性 的 平均 提高 率 


平均 提高 


参数 


对 比 PF 对 比 FA-PF 


Q-1 22.6% 7% 


Q=0.1 5.7% 3.8% 


表 3 数据 可 知 ，Q=1 时 ， 本 文 算法 准确 率 较 最 新 的 FA- 
PF 算法 提升 7%， 较 PF 算法 有 22.6% 的 提升 ， Q=0.1 时 ， 本 文 
算法 的 准确 性 较 FA-PF 和 PF 有 3.8% 和 5.7% 的 提高 。 这 是 由 于 
本 文 算法 为 确保 滤波 器 的 精度 ， 在 萤火虫 的 迭代 寻 优 过 程 中 引 
入 了 最 新 的 观测 值 信息 。 

2) 算法 跟踪 效果 测试 
在 相同 过 程 噪声 , 不 同 粒 子 数目 N=50、N=100 条 件 下 运行 
本 文 算法 与 PF、FA-PF 算法 并 记录 三 种 不 同 算法 跟踪 的 状态 估 
计 值 及 误差 绝对 值 数据 ， 结 果 如 图 4、5 所 示 。 

a) 当 N-50, Q-1 时 ,实验 的 估计 值 及 误差 绝对 值 分 别 如 图 
4(a) (b) 所 示 。 


| Fest 


,| |*PF 估 计 值 | M WU 
:FA-PF 估 计 值 | 3 i Pow 
* 本 文 算法 估计 值 
"i : : i i i a 
时 间 步 长 


(a) 滤 波 估 计 值 


时 间 步 长 
(b) 误 差 绝 对 值 


吉 果 与 其 他 算法 对 比 


图 4 N=50 时 本 文 算法 的 跟踪 结 
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王 航 星 ， 等 : om 半径 的 营 火 法 的 粒子 滤波 


b) 当 N=100，Q=1 时 ， 实 验 的 估计 值 及 误差 绝对 值 分 别 如 


图 5(a) (b) 所 示 。 
f 
TI 
A | 
20 
LA J 
ji f H 
Akt | | 
“PF 估计 什 i| xy EV 
“FA-PF 估 计 值 yoo 6 DON 
。 本 文 算法 估计 什 S 
IE. ^ MEER 
(a) 滤 波 估 计 值 
“PF 误差 
+ —FAPFRJX 
* 本 文 算法 误差 


长 
(b) 误 差 绝 对 值 

图 5 N=100 时 本 文 算法 的 跟踪 绰 
从 图 4(a)、5(a) 可 以 看 出 , 本 文 算法 在 不 同 的 粒子 数目 下 的 


吉 果 与 其 他 算法 对 比 


滤波 估计 值 要 比 PF 及 FA-PF 算法 更 加 接近 于 真实 值 ， 反 映 在 
图 4(b)、5(b) 中 即 本 文 算法 的 误差 绝对 值 要 小 于 PF 及 FA-PF 算 
法 ， 本 文 算法 的 准确 性 在 粒子 数目 不 同 的 情况 下 均 高 于 PF 及 
FA-PF 算法 。 
这 主要 是 由 于 本 文 算法 在 萤火虫 算法 中 采用 了 融合 粒子 观 
则 值 信 息 的 相对 亮度 计算 方式 ， 在 粒子 集 的 迭代 寻 优 过 程 中 根 
据 相对 亮度 更 加 准确 地 调整 粒子 的 位 置信 息 ， 这 使 得 粒子 更 加 
准确 地 分 布 在 其 真实 值 附 近 ， 从 而 算法 的 准确 性 得 到 了 一 定 的 
提高 。 
4.2.2 实时 性 测试 
在 软 硬 件 条 件 相同 的 条 件 下 ， 在 不 同 粒子 数目 、 不 同 过 程 
噪声 条 件 下 运行 不 同 算 法 并 记录 算法 运行 时 间 进 行 对 比 ， 结 果 


ue 


对 比 见 表 4。 
表 4 不 同 粒子 数目 下 几 种 算法 运行 时 间 对 比 
参数 PF ”FA-PF 本 文 算法 对 比 FA-PF 缩减 率 
N=50 0.0101 0.0264 0.0227 14.02% 
Q-1 N=100 0.0175 0.0629 0.0585 7.0096 
N-200 0.0404 0.1685 0.1515 10.09% 
N=50 0.0093 0.0194 0.0177 8.7% 
Q=0.1 N=100 0.0152 0.0457 0.0431 5.7% 
N-200 0.0323 0.01241 0.1241 8.8% 
从 表 4 可 以 看 出 ,在 不 同 粒子 数目 \ 不 同 过 程 噪声 条 件 下 ， 


本 文 算法 的 运行 时 间 较 PF 算法 存在 一 些 延 长 ， 但 结合 表 2 数 
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据 ， 本 文 算法 的 准确 性 较 PF 算法 有 大 幅度 的 提高 。 这 是 由 于 
本 文 算法 在 重 采样 前 使 用 了 粒子 集 的 迭代 寻 优 操作 ， 使 多 数 粒 
子 的 权 值 都 较 高 ， 所 以 本 文 以 运行 时 间 加 长 为 代价 换 来 准确 性 
的 提高 。 并 且 在 不 同 过 程 噪声 条 件 下 ， 三 种 算法 的 运行 时 间 随 
着 粒子 数目 的 增加 都 有 一 定 的 延长 。 

通过 表 4 数据 计算 可 得 在 不 同 的 过 程 噪声 下 ， 本 文 算法 的 
运行 时 间 相 对 于 FA-PF 的 平均 缩减 率 如 表 5 所 示 。 
RS 不 同 过 程 噪声 下 ， 本 文 算法 运行 时 间 的 平均 缩减 率 


参数 平均 缩减 率 
Q-1 10.494 
Q-0.1 7.7% 


dé 4 和 5 中 数据 可 以 看 出 ， 在 不 同 粒 子 数目 以 及 不 同 过 
程 噪声 的 情况 下 ， 本 文 算法 的 运行 时 间 较 FA-PF 算法 都 存在 一 
定 程度 的 缩减 。 在 Q-1 及 Q-0.1 时 ， 本文 算法 的 运行 时 间 对 比 
FA-PF 算法 有 10.4% 和 7.7% 的 平均 缩减 。 

这 是 由 于 在 使 用 萤火虫 算法 进行 粒子 间 的 相互 吸引 完成 迭 
代 寻 优 过 程 中 ， 本 文 提 出 并 使 用 了 自 适 应 吸引 半径 来 控制 粒子 
的 吸引 范围 ， 从 而 大 大 减少 了 参与 计算 的 粒子 数目 ， 使 得 算法 
的 运行 时 间 有 了 进一步 的 缩减 。 
43 定性 分 析 

为 了 直观 地 体现 出 本 文 算法 的 跟踪 性 能 ， 在 主流 的 视频 肌 
Visual Tracking BenchmarkP20 上 选取 视频 进行 实验 。 在 跟踪 过 程 
中 ， 采 用 颜色 直方 图 作为 跟踪 特征 ， 取 粒子 数 N=100。 将 本 文 
算法 与 PF 算法 的 跟踪 效果 进行 对 比 ， 并 对 跟踪 结果 的 部 分 帧 
进行 截图 ， 如 图 6~9 所 示 。 
在 Coke 视频 中 ， 跟 踪 场 景 较为 简单 ， 跟 踪 目 标的 形态 也 
未 发 生 大 幅度 变化 ， 跟 踪 目 标 作 无 规律 的 圆周 运动 ， 是 典型 
的 非 线 性 系统 下 的 目标 跟踪 。 图 6 为 视频 序列 的 完整 场景 图 ， 
图 7 为 PF 算法 跟踪 视频 序列 结果 , 图 8 为 文献 [14] 中 提出 的 的 
FA-PF 算法 跟踪 视频 序列 结果 ， 图 9 为 本 文 算法 跟踪 视频 序列 
结果 。 
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图 6 Coke 视频 完整 场景 图 
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图 9 本 文 算法 跟踪 Coke 视频 
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录用 定稿 
图 7-9 可 明显 地 看 出 , PF 算法 的 粒子 分 布 很 散乱 ,粒子 
分 布 不 均匀 使 得 对 跟踪 有 意义 的 粒子 数目 大 大 减少 ， 这 是 导致 


PF 


或 ， 


wu 


附近 


法 更 为 集中 ， 本 文 算法 的 粒子 分 布 最 为 
边界 框 的 中 心 位 置 ， 这 是 由 于 本 文 算法 利用 改进 的 萤火虫 算法 
对 粒子 集 进行 了 迭代 寻 优 操作 。 


在 第 80 帧 和 第 131 帧 ， 由 于 散乱 的 粒子 无 法 集中 到 高 似 然 区 


跟踪 。 


算法 跟踪 失败 的 重要 原因 。FA-PF 算法 的 粒子 分 布 较 PF 算 
中 ， 主 要 集中 在 跟踪 


Yit 


在 第 38 帧 ， 跟 踪 目 标 运动 到 场景 的 后 方 发 生 了 部 分 遮挡 ; 


这 些 原因 都 导致 PF 算法 跟踪 失败 。 相 比 于 基本 实现 准确 


REXHJ FA-PF 算法 而 言 ， 本 文 算法 使 用 改进 的 萤火虫 算法 有 效 
也 对 全 部 粒子 进行 迭代 寻 优 ， 使 粒子 集 尽 可 能 多 地 分 布 在 目标 


Tp 


， 表 示 跟 踪 有 贡献 的 粒子 数量 较 大 ， 从 而 实现 更 加 准确 的 


除 此 之 外 ， 本 文 算法 由 于 引入 了 自 适 应 吸引 半径 ， 有 效 控 


制 了 参与 寻 优 的 粒子 数目 ， 缩 减 了 迭代 寻 优 过 程 。 由 表 4 数据 


可 知 ， 算 法 运行 的 时 间 也 有 了 7% 的 缩减 。 
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粒子 滤波 算法 需要 大 量 的 粒子 才能 保证 算法 的 准确 性 ， 随 
值 的 不 断 削减 ， 用 于 跟踪 的 粒子 数目 越 来 越 少 ， 需 要 对 粒 
进行 重新 采样 。 但 这 会 严重 削减 粒子 的 多 样 性 ， 甚 至 导致 


法 对 粒子 滤波 算法 进行 优化 。 


结合 粒子 的 观测 值 信息 使 粒子 集 分 布 在 全 局 最 优 的 区 域 ， 保 证 


只 剩 下 一 种 粒子 ， 导 致 跟踪 失败 。 本 文 使 用 改进 的 萤火虫 


根据 两 种 算法 的 特点 ， 对 萤火虫 的 亮度 公式 进行 了 改进 ， 


了 算法 的 准确 性 ， 同 时 提出 了 自 适 应 吸引 半径 参数 ， 避 免 了 粒 
子 数目 增加 时 计算 量 的 大 幅 增长 ， 降 低 了 算法 的 运算 复杂 度 。 
定量 分 析 和 定性 分 析 均 表明 ， 本 文 算法 有 效 地 减少 了 跟踪 过 程 


中 月 


确 的 跟踪 ， 同 时 本 文 算法 的 运行 时 间 也 相对 于 FAPF 有 


f 需 的 粒子 数目 ， 在 粒子 数目 较 少 的 条 件 下 依然 能 够 实 
J 


的 缩减 。 
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